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Abstract 

We study chiral symmetry breaking in QED when a uniform external magnetic field is present. 
We calculate higher order corrections to the dynamically generated fermion mass and find them 
to be small. In so doing we correct an error in the literature regarding the matrix structure of the 
fermion self-energy. 







In Ref. the effects of an external magnetic field on chiral symmetry breaking was stud- 
ied using quantum electrodynamics as the model gauge field theory. The Schwinger- Dyson 
equation for the fermion self-energy in the quenched, ladder approximation was expressed in 
terms of Ritus' representation [2] for the exact fermion propagator in a constant magnetic 
field and an approximate solution for the dynamically generated fermion mass was found. 
The infrared fermion mass found in Ref. [l| was consistent with that obtained from other 
approaches ^. 

In the present work we will examine how robust this approximate result is by carrying out 
the calculation to the next leading order of the approximation. In the course of our study, 
we discover an inconsistency in Ref. jl| regarding the form of the dynamical fermion mass 
matrix. We obtain the correct matrix structure and show that, despite this inconsistency, 
the infrared dynamical fermion mass found in Ref. |l| remains correct. 

Let us begin with the Schwinger-Dyson equation for the fermion self-energy in the 
quenched, ladder approximation (see Eq. (34) in Ref. jl|): 
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We adopt the notation and convention of Ref. [1| and will not explain them here. Since 
we are interested in finding a solution for the infrared fermion mass, the calculation will 
be simplified by considering the p — limit in Eq. (^. Let us remind the reader that the 



4-momentum p has the components {po, 0, —sgn{eH) ^^J 2\eH \ k, p^) . Consider first the k ^ 
limit. Due to the presence of 6kk' on the left hand side of Eq. (^, letting k = implies k' = 
which in turn requires n = n' = and a = a' = sgn{eH). This is because the quantum 
numbers n, k, and a are related by the relation 

n = n{k, a) = k + ^sgn(eif) - ^ (2) 

with the allowed values n = 0, 1,2, k = 0, 1,2, and cr = ±1, and similarly for n' = 
n{k',a'). In this limit, the J polynomials in Eq. (P) become 

Jon"{qL) = [i sgn{eH)q_L\''" , 

Jn"o{q±) = sgn{eH)q,_f'. (3) 



and Eq. (jT)) now reads 
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where A is now understood to be A = A(sgn(eif)), p\\ = {po,P3), and the 4-momentum p" 



in the integrand is given by p" = (po — Qo, 0, —sgn{eH) y 2\eH\k" , p^ — q^). 

By transforming to polar coordinates for the integrals over qi and q2'. J dqidq2 = 
I qLdq±dip, the integration over the angle ip can be carried out to yield 

J rf<^e'«^-(-^^)(""-"")'^ = 2n5n"n" = 2Ti5-„n,„. (5) 

This allows us to rewrite Eq. as 

. .^p|e^|)£i:-L/||,e-»l(z|)!: 
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where A" = A(a")- 

The fermion self-energy is expected to have the form = Z{p)'y ■ p + where 

^a{P) is a matrix representing the dynamically generated fermion mass. An ansatz was 
made in Ref. jl| that was proportional to the unit matrix. An approximate solution 
for was then obtained by keeping only the dominant k" = term on the right hand 
side of Eq. ©. In this case, n" = and a" = sgn(eif). However, in the calculation of 

n ^ 

Ref. PJ , the spin summations in Eq. (|T]) were carried out before the infrared limit was taken 
(see Eqs.(44)-(47) there), and contributions from both a" = sgn(eiJ) and a" = — sgn(ei7) 
were included for the k" = term, thus leading to the incorrect conclusion that was 
proportional to the unit matrix. We shall see below that should be proportional to 
the matrix A(sgn(e-ff)), which is equal to diag{l, 0,1,0) for sgn(eiJ) = +1 and equal to 
diag{0, 1, 0, 1) for sgn{eH) = -1. 

If we let ^a{p) = rn{p)A, where m{p) denotes the dynamically generated fermion mass, 
and keep only the k" = contributions, Eq. © becomes 



Z(p||)7-P||+m(p||)A ^ ^^'(2|ei/|)/^^L.-(l-0 ^ 
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where pj| = p\\—q\\- Because 7 -pjl commutes with A, (— [l + Z]7-pj| +mA)([l+Z]7-pj| +mA) = 
[1 + + m^A, which is a diagonal matrix that also commutes with 7|j' = (7°, 7^). We 

can therefore rewrite Eq. (jZj) as 
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If we work in the Feynman gauge = 1) and use the property that 

A7'^A = A7ff A = 7ffA = A7ff, (9) 



we find 
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[l + Z(pj[)]2(pj[)2 + m2(pj[)A 
— 2m(pj|) 
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and Eq. (jHJ becomes 



Z(p„h.p„+™(p,)A== -'e'Pkg|)/(^-^ |,^^(^pp(p,,;.^„.(p,,) A. (11) 

It follows that Z{p\\) vanishes in the Feynman gauge, as was found in Ref. and one 
obtains the Schwinger-Dyson equation for the dynamical fermion mass: 

Note that both sides of this equation are proportional to the A matrix, which justifies our 



ansatz for the matrix structure o: 
Following the authors of Ref. 



l|, we seek a solution for m in the infrared limit {p\\ — > 0) 
and approximate wi(pj|) that appears in the integrand on the right hand side of Eq. (fT^ by 
its infrared value. This leads to the gap equation 
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where we have made a Wick rotation to Euchdean space. This is precisely the gap equation 
(52) in Ref. Q], which has a solution of the form 



mo ~ a^J\eH\e (14) 

where a and h are real positive constants of order one and a is the fine structure constant. 
The subscript indicates that this is the lowest order approximation for m. As noted in 
Ref. Ill, this solution is applicable for small a. 

Our result confirms the finding for the lowest order solution to the Schwinger- Dyson 
equation obtained in Ref. We also obtain the correct matrix structure for S^, which is 
proportional to the matrix A(sgn(eiJ)). 

Let us consider now the higher order corrections for the dynamical fermion mass m. In 
the Feynman gauge, Z{p\\) = and the Schwinger-Dyson equation Q may be written in 
the infrared limit as 
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where we have set S^(0) = mA and separated the k" = term which contributes to the 
lowest order solution rriQ. We shall estimate the next order corrections to m by replacing 
m{p") in the integrand with the lowest order solution mo. We need therefore to evaluate 

1 / i-qir" 
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Aj'^A"^^ ^A"7^A, (16) 



where p" = {—qo, 0, — sgn(eif )-y/ 2|eif|/c", — gs). 

We shall consider the corrections to m coming from the following contributions. First, 
for n" = 0, there is a subdominant k" = 1 (and hence a" = —sgn{eH)) term which was not 
included in mo. Second, there are two terms for n" = 1: one for k" = 1 and a" = sgn(eif), 
the other for k" = 2 and cr" = —sgn{eH) . 

For k" 7^ 0, (7 ■ p") no longer commutes with A, and the expression for (7 ■ p" + moA)~^ 
is no longer so simple and depends on sgn(eif). We shall present our calculation below 
for the case sgn(eif) = +1. However, the final result is applicable for either sign. For 
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sgn(eiJ) = +1, one finds that [7 ■ p" + moA(l)] ^ is given by 
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where k = \j2\eH\k" , X = qo + q^, and P = qo — q^. Fortunately, this simphfies a great deal 
after we compute A7^A"(7 ■ p" + moA)~^A"7^A. 
Since a" = ±1, and 

A(l)7iA(l)=0 = A(l)7fA(-l), (17) 
where 7;^ = (7-^,7^), we need only evaluate 
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and 
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where (T3 denotes the Pauli matrix. The off-diagonal elements in Eq. might seem to 
be a problem. However, more careful inspection shows that they are odd in go and q^ and 
will vanish upon integration over these variables, thus reducing the matrix to a diagonal 
matrix proportional to A(l). It should be stressed that both matrix structures appearing in 
Eqs. (|18p and (jl9j) are proportional to A(l), consistent with the left-hand-side of Eq. 
This shows that is proportional to the matrix A{sgn{eH)) for all higher order terms in 

Eq. (uni). 

Putting all the pieces together, using Eq. (fTHj) for the k" = 1, a" = 1 term and Eq. (fT^ 
for the other two terms corresponding to a" = —1 and k" = 1,2, we find the next order 



corrections for the dynamical fermion mass to be 



(m-mo)A ~ ze^(2|eif|) 
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In other words, 

m — mo 
mo 



TT 



gj + gi I2|eff|(l + gj)2 + m2g2 
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_2|eif|(l + gJ)2 + m2g2 2|eiJ| (2 + gf )2 + mggj 

where we have performed a Wick rotation to Euchdean space and carried out the integration 
over the polar angles on the q\\ plane as well as on the q± plane. 

We see that the fractional correction to m is quite small, being proportional to (a/vr). If 
we substitute the expression for mo found in Eq. (HH), with a and h set to be one, we can 
evaluate the integrals numerically to obtain the estimate 

m — mo 



mo 



a 



TC 



(0.4275 - 0.2161 - 0.1628) 



a 

0.0486 ( - 

TT 



(22) 



a very small number indeed. 

In summary, we have reexamined the Schwinger- Dyson equation analysis of chiral sym- 
metry breaking in QED in the presence of a uniform external magnetic field, confirming the 
lowest order result for the dynamical fermion mass obtained in Ref. ^l? and finding the cor- 
rect matrix structure for the fermion self-energy. We have also calculated the higher order 
corrections to the dynamical fermion mass and found them to be small, thus validating the 
approximation scheme used in the calculation. 
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